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Titration of methacrylic acid / ethyl-acrylate copolymers is studied experimentally and theoret- 
ically. At low salt concentrations, this polyacid exhibits a plateau in the titration curve below the 
neutralization point. The plateau has been often attributed to a first-order phase transition associ- 
ated with polymer conformational changes. We argue that the specific shape of titration curves of 
hydrophobic polyelectrolytes is due to electrostatics and does not necessarily require a conformation 
change of the polyelectrolyte chains. We calculate the free energy at the mean-field level and its 
first-order (one loop) correction using a loop expansion. The latter is dominated by Debye-Hiickel- 
like charge-charge correlations as well as by correlations between dissociation sites along the polymer 
chain. We show that the one-loop corrections to the free energy lead to titration curves that agree 
with experiments. In particular, the model explains the decrease of the pH at the plateau when the 
polymer concentration is increased or when salt is added to the solution. 



I. INTRODUCTION 

In recent years the search for environment-friendly materials has promoted the development of numerous water- 
soluble polymer applications. Polymers can be made water soluble by making them compatible with the strong polar 
environment of the aqueous media, e.g., by introducing either charges or strong dipoles on the chains [[l] ||. 

Most charged polymers (polyelectrolytes) have a hydrophobic backbone. This hydrophobicity induces an effective 
attraction between monomers which competes with the Coulomb repulsion between charges. As a result, hydrophobic 
polyelectrolytes exhibit complex behavior including conformation changes, macro- and meso-phase separation, self- 
association and aggregation 

In spite of considerable theoretical and experimental effort, many questions remain open. In particular, the role of 
the coupling between hydrophobic attractions and long range Coulomb interactions on the physico-chemical properties 
of polyelectrolyte solutions is not fully understood. In this paper, we consider the role of such coupling in solutions 
of weak polyacids. 

Weak acid monomers (denoted HA) can undergo dissociation of the type 

HA # H+ + A" 

where A - is the charged monomer attached to the backbone, while the dissociated H + charge dissolves into the 
solution. The dissociation/association is an equilibrium process satisfying detailed balance. The ionization degree 
determines the effective amount of charge on the polyelectrolyte chain, and it depends on the pH of the solution. 
At low pH the polymer is weakly charged while at high pH a larger fraction of monomers is dissociated and the 
polymer charge saturates to its maximal value. The most visible consequence is the solubility in water: hydrophobic 
polyacids can become water soluble at high enough pH, where the polymer charge is strong enough to overcome the 
chain hydrophobicity. In contrast to low molecular weight acids, the charged groups of polyacids are correlated as 
they are linked together along the chain. Indeed, the dissociation of one acid group is correlated in a complex way to 
the position and number of other charged groups on the chain. As a result, when the amount of charge on the chain 
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varies, the chain conformation is affected, and in turn, influences the dissociation of other groups yg0,§. Studies of 
the interplay of these competing effects has attracted a large amount of experimental |p|-[ll[| and theoretical |l2]-|l4|] 
interest. Additional motivation for these efforts is related to the use of water-soluble polyelectrolytes and in particular 
alkali-swellable polymers in many industrial applications, such as coatings, food and cosmetic industries (TsIJk] ]. 

An elementary and wide spread experimental tool to characterize polyacids is to perform titration experiments. In 
these experiments, a strong base like NaOH is added to a solution of weakly charged polyelectrolytes. The pH of the 
solution and the equilibrium dissociation depends not only the amount of added base (like for ordinary acids 
but also on the polyelectrolyte concentration and presence of salt |l7|Jl^ ]. The overall shape of the titration curve 
depends on the nature of the polyelectrolyte which is titrated ]l7| ]. In the case of poly (acrylic) acid, the pH increases 
steadily with the dissociation degree. In contrast, the titration curve of poly(methacrylic) acid shows a maximum at 
low dissociation. Such a non-monotonous dependence has been related to a conformational transition. The case of 
hydrophobic polyelectrolytes is still more intriguing: the titration curve shows a large plateau where the pH is almost 
constant before the neutralization point is obtained. It has been argued in the literature that this plateau could be 
associated with a first-order phase transition between collapsed and swollen states of the polyacid chains j?],^0| . 

In this paper we show that this peculiar shape of the titration curves can be explained without the need to rely 
on conformational phase transition of the polyacid chains. By including the effect of correlations of charges along 
the chains in the free energy, titration curves are calculated and show a behavior similar to the curves obtained in 
experiments. In particular, we explain how the pH depends on the polymer concentration and the amount of added 
salt. This dependence is special to polyacid solutions and is much weaker in monomeric acid solutions. 

The paper is structured in the following way: in the next section the mean-field free energy, its one-loop correction 
and the resulting titration equations and curves are presented. Then, the experimental measurements of the titration 
of methacrylic acid / ethyl-acrylate copolymers (MAA-EA) are discussed in Sec. Ill, and the comparison between 
theory and experiment is presented in Sec. IV. The full formalism relying on a field theoretical approach will be 
presented in a forthcoming publication. 



II. THEORY 



We present first the free energy leading to titration curves of weak polyacids in presence of added salt. Our system 
consists of four dissociating species: water, methacrylic acid monomers (denoted HA), NaOH titrating base and NaCl 
salt. The dissociation reactions are written as: 

H 2 O^H+ + OH-, 
HA ^ H+ + A", 
NaOH — > Na + + OH~, 
NaCl->Na+ + Cl-. 

The partial dissociation of water and acid monomers is accompanied by an energy cost of breaking the molecular 
bond. They are denoted as Ai and A2, respectively (in units of the thermal energy k^T), and they are related to the 
mass action law as will be detailed below. Since NaOH is a strong base it is fully dissociated. Similarly, the salt 
is fully dissociated. 

The polymer used in the experiment is a statistical copolymer composed randomly of two monomers: methacrylic 
acid and ethyl-acrylate. The total polymerization index is denoted as N. A fraction / = 1/3 of the monomers are 
composed of the methacrylic acid, i.e. they are ionizable. Since the charged monomers are distributed uniformly 
along the chain, we assign a partial charge fe to each monomer ( "smearing" the charges along the chain) . Note that 
fe is the nominal charge of each monomer, while the actual charge is related to the partial dissociation of the acid 
monomers. 

The concentration of monomers in the solution is denoted c m while the concentrations of the small ions are denoted 

CH) COH, CNa, Cci- 

The base concentration added to the solution is denoted cb and its ratio to the MAA monomer concentration is 
defined as the degree of neutralization 7 

where 7 varies between zero (no added base) and infinity (large base excess). It is also useful to define the degree of 
ionization a of the acid 
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„ = [A'] = [At] 

[HA] + [A-] fc m 1 ' 

where [A - ] and [HA] are the concentrations of dissociated and non-dissociated monomers, respectively. This degree 
of ionization is related to the pH=— log 10 ch of the solution via 

pH = pK A + log 10 —— (3) 
1 — a 

In addition, due to charge neutrality, a and 7 are related via 

[H+] - [OH-] 

a = 7 + J 1 (4) 

/Cm 

The conservation of mass implies the following relation between the ion concentrations: 



CNa = C B + CCJ 

Note that the Cl~ ions come only from the salt, while the Na + come both from the base and the salt. Similarly, the 
OH - ions come from the dissociation of water and base, whereas the H + ions come from the dissociation of water 
and monomeric acid. 



A. Free energy 

We denote the full free energy as F = Fq + AF, where Fq is the mean field term and AF is the one-loop correction. 
Below we discuss the two terms of the free energy separately. 



1. Mean- field free energy 

In a separate work, the mean-field free energy (per unit volume and in units of fceT) is shown to be 

(3F = ^(logC m W m -l) + Y C J ( lo S C W J - !) - | C m + 7T C m ~ A C m 

j=H,OH,Na,Cl " ' 

- (coh - Cb) Aj - (c H + c B - c OH ) A 2 + (ch + c B - c OH ) /3e</? - /c m log(l + e^°) (5) 

It turns out to be identical to the usual Flory-Huggins theory of polymer-solvent mixtures. The factor 1/N in the 
first term accounts for the reduction of translational entropy of chains of N monomers. The second term represents 
the entropy of mixing of the small ions. The parameter Wj is the molar volume of the j species. The next two terms 
represent the short range monomer-monomer interactions. The hydrophobic effect is modeled by an attractive second 
virial coefficient v (having dimension of volume), and a repulsive third virial coefficient w (having dimension equal to 
volume squared). It is introduced to avoid collapse of the chain. The next term represents the chain conformational 
entropy, where Ao is the conformational entropy per monomer of a free Brownian chain. 

The following two terms represent the energy cost of dissociation of water molecules and acid monomers. We define 
by Ai the energy loss for each dissociation of a water molecule 



H 2 ^ H+ + OH~ 

and by A2 the energy loss for each dissociation of an acid molecule 

HA ^ H+ + A" 

The number of OH~ ions coming from dissociation of water molecules is equal to the difference between the total 
number of OH~ ions and those coming from the NaOH base. Similarly, the number of dissociated acid groups is equal 
to the difference between the total number of H + ions and the number of H + ions coming from the water. 

Finally, the last two terms account for the electrostatic energy, ipo being the electrostatic potential in the solution. 
The first one is the electrostatic energy of the small ions, while the second is the electrostatic free energy of partially 
dissociated monomers Ell] . 
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2. One-loop correction to the free energy 



The one-loop correction to the mean-field free energy AF will be presented in detail in a forthcoming publication. 
It is obtained by integrating over the quadratic fluctuations of the concentration fields. The correction term to the 
free energy is given (up to an additive constant) by 

E(<z) = [1 + (wc m - v)c m ND(r,)] (q 2 + n 2 ) + A7rl B f 2 A 2 ND(r ] )c m (7) 
V = ^a 2 q 2 N (8) 

where the Debye- Hiickel screening length k _1 depends on the total concentration of small ions ci and the effective 
charge of the chain fA(l — A)c m 

k 2 = 4tt/b(ci + fA(l - A)c m ) = 4nl B c cS (9) 
c\ = c H + c h + c B + 2c C i (10) 

Ib = e 2 jek^T ~ 7 A is the Bjerrum length. The dissociation fraction of the acid monomers is given by 

(11) 



1 + e^ e v° 

The Debye function D(r/) ||] entering eq. (0) is given by 

2 e~ v - 1 

D( V ) = -(1 + ) (12) 

rj Tj 

In addition to the small ion contribution to the electrostatic screening, eq. @ includes a term proportional to 
A(l — A), where A is defined in eq. ( pr] ) above. This term accounts for changes in the dissociation degree of monomers 
depending on the local electrostatic potential. At small A values, A(l — A) w A, and the polymer contribution is the 
same as that of disconnected monomers. However, as A — > 1 there is a substantial reduction (by the factor 1 — A) of 
the effective polymer charge that contributes to screening J§|,^J . This reduction can be understood in terms of charge 
correlation along the chain. Dissociated monomers prevent further dissociation of other monomers. 

In our model, the concentrations CN a , cci and c m are fixed by the amount of NaCl salt, NaOH base and polymer in 
the solution while the concentrations of dissociated H + and OH~ ions depend on the degree of dissociation of water 
molecules and acid monomers (through the electrostatic potential i^o)- Thus ch, coh and (p$ are variational parameters, 
determined by the requirement that the free energy F = F a + A_F is an extremum. In titration experiments ch can 
be monitored directly through the pH of the solution. 



3. Simplified free energy 

In order to better understand the main contributions to the correction AF let us return to the expression for S((j), 
eq. (0). First, we note that the Debye function is bound between and 1 (see also eq. (^7|)). For typical values of the 
physical parameters (see also Sec. IV), 

v ~ 25A 3 w ~ 100A 6 

N ~ 10 4 c m ~ ImM (13) 
it is easily seen that for long chains (N ^> 1), the Debye function can be approximated by ND(r]) ~ 12/a 2 q 2 , and 

\wc m - v\ Cin ND{r)) « 1 (14) 
Thus, the effect of the solvent is small and will be neglected in this section 

A ^"^/0 1Og ^ + K2 + Z ] (15) 
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where 

z = 4i:l B f 2 A 2 ND(7 1 )c lri (16) 
For the same range of parameters k 2 3> z and AF can be expanded to first order in z, AF = AFq + AF±: 



* F ° = ~o T^M^ + K 2 ] (17) 



2 7 (2tt)3 

The integral can be calculated by including a cut-off at large q, which cancels out the Coulomb self-energy, yielding 

AF = - (18) 

This result is the well known Debye-Hiickel correlation energy of an electrolyte, because the polymer contribution 
is negligible and does not appear in this leading term. The correction AFi is given by 

AF, ~ 1 / Z ~ 6fM 2 ^^c (19) 

This correction is due to correlations along the chain between dissociated monomers as can be seen from the 
following simple argument. Consider the screened electrostatic interaction between charged monomers on a single 
infinite chain. For a specific chain configuration the Coulomb energy per monomer is 

f2 /|2„2 



where rj is the spatial distance between the j = monomer and another j = ±1,±2, ... monomer. For a Gaussian 
random walk, the typical (most probable) distance between monomers is given by 



fa) = [jj a (21) 

Assuming that rj can be approximated by (rj) in eq. (|20|), and na <C 1, the sum over j can be replaced by a continuous 
integral leading to eq. (|l9|). 

B. Titration equations 

The free energy depends on the species concentrations: coh, CNa, ch, c m , cci. However, the concentration of 
monomers, Na + and Cl~, is fixed by the amount of polymer, base and salt added to the solution. On the other 
hand, the concentration of H + and OH~ ions is determined self-consistently by minimizing the full one-loop free 
energy F = Fq + AF, eqs. (jfj), (g), with respect to ch,coh and ipo. The resulting equations of state determine the 
dependence of the pH of the solution on the other system parameters. 

log(cH^ H ) + petpo - (3A 2 + 2nl B h = (22) 
log(coH^OH) - /3e<A, - (3A X + f3A 2 + 2nl B h = (23) 

ch + cb- coh - fcnA^l - 2irl B {l - A)({1 - 2A)h + 2ANfI 2 )^j = (24) 

where the quantities I\ and I 2 are defined as: 

f A k 2 dk l + Nc m (wc m -v)D(r ] ) 

J l = / ~~n 2 ( 25 ) 



Jo 2tt 2 E(fc) 

, 2 =f^m (26) 



o 2tt 2 S(fc) 



and A is a short distance (large q) cut-off. 
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With the full expression (12) of the Debye function D(rf), the integrals in (25), (26) cannot be performed analytically. 
However, for the values of the parameters used in the experiments, namely polymerization index N ~ 10 — 10 5 and 
monomer length a ~ 5A, it is possible to use a simplified form for the Debye function: 

D(n) ~ ^—r (27) 

K " 1 + rj/2 v ; 

As can easily be seen by comparing eq. ( pTj ) with the exact expression (p^), the above form has the right behavior 
both at small and large values of rj. Within this approximation, the integrals are given by: 

7 1 R 2 k 2 k + k_/2 + (\ + 1)k 2 + u 

2nR 2 (k + + k_) 

h = 2nR 2 (K + + K _) (29) 



where we have used the notation: 



A = Nc m (wc m - v) (30) 
u = Anl B A 2 Nf 2 c m (31) 
R 2 = Na 2 /6 (32) 



and 



4 - jp ( A + 1 + ± \A A + 1 ~ r2k2 / 2 ) 2 - 2uR2 ) ( 33 ) 

Note that in eq. ( pq ) we have omitted a term, equal to A/27T 2 , which exactly cancels the Coulomb self-energy. 
With these definitions, eqs. ( p2f ) and ( |23| ) can be recast in the form: 

c H -coh - 10- 14 e- 4TdB/l (34) 

which shows the change induced by fluctuations on the water dissociation constant and 

c H = ^—t^ 10~ pKA c - 27Tlsh (35) 
J\ 

These two equations, together with eq. (p4|), which expresses charge neutrality at the one- loop level, are solved 
iteratively. The numerical solution is obtained by using the mean field values as a starting point for the iterations, 
and convergence is usually achieved after a few iterations. 



C. Structure Function — S(q) 

In scattering experiments the structure function S(q) is readily obtained. It is related to the Fourier transform of the 
various density-density correlations in the sample. For example, we can regard the monomer-monomer correlations: 

S(q) = (Sc m (q)Sc m (q))/c m (36) 

Since the one- loop expansion takes into account Gaussian fluctuations, it can also be used to calculate the structure 
function. The result is 

„-i n 1 2 f 2 A 2 4irl B c m 

S {q) = NW) + m ~ m+ <Z 2 + * 2 (3?) 
The inverse structure function S~ 1 (q) is the energy penalty associated with density fluctuations at a wavenumber q. 
A minimum in S~ 1 (q) corresponds to the strongest fluctuating wavenumber q* . As a result, incoming radiation at this 
wavenumber interacts most strongly with the sample and a peak appears in the structure function and consequently 
in the scattering intensity. 

An instability appears when S~ 1 (q) vanishes, corresponding to a divergence of the peak in the structure function. 
When the instability appears at a finite q = q* the system undergoes a meso-phase separation and becomes spatially 
modulated. If, however, the instability is at q = the system undergoes a macrophase separation. In our case this 
occurs when 

1 i 2 A 2 r 

uc m = wc m + — H (38) 

N c off 

Recall that c c s and A depend on the pH and the degree of neutralization. 
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III. EXPERIMENT 



A. Preparation of samples 

The copolymer polyelectrolyte chains used in this study are prepared by standard emulsion polymerization tech- 
niques using neutral ethyl-acrylate (EA) monomers and methacrylic acid (MAA) monomers. The methacrylic acid is 
a weak acid with pKa=4.5. The weight fraction of MAA in the copolymer is equal to 0.35 taken to be / = 1/3 in 
the theoretical section. The emulsion polymerization is performed under starved monomer conditions. Under these 
conditions, the MAA and EA monomers are evenly distributed along the polymer chain. From size exclusion chro- 
matography, we have estimated the molecular weight of the chains to be of the order of 106 daltons, corresponding 
to a polymerization index N ~ 10 4 . The MAA-EA copolymers precipitate from the aqueous solution, as they are 
insoluble in water. The precipitate is carefully washed by ultrafiltration in order to remove surfactants, unreacted 
monomers and initiators. This cleaning procedure is stopped when the resistivity of the water flushed through the 
separation membrane is that of pure water (18.2 Mfi/m). The solid content of the stock solution is then determined 
accurately by drying and weighing. 

B. Titration measurements 

Polymer solutions at different weight concentrations are prepared by mixing weighted amounts of the stock solutions 
with de-ionized water. The dissolution of CO2 is prevented by carefully de-gassing the solution with nitrogen. Each 
of the polyelectrolyte solutions is neutralized by an NaOH solution with a molar concentration ranging from 0.1M to 
2M depending on the polymer concentration. Upon neutralization, the methacrylic acid is neutralized and repulsive 
forces due to the negative charges cause the chain to expand, resulting in the progressive solubilization of the polymer 
chains. As a result the solution becomes transparent and its viscosity increases. In the following, we shall characterize 
the neutralization of the polymer chains by the degree of neutralization, 7 = CB/(/c m ), which is the ratio of the 
amount of added base to the amount of available acids groups. Titration experiments are performed using a pH-meter 
(Metrohm 691) with a combined glass electrode. The measurements are made at 20°C with constant stirring under 
a nitrogen atmosphere. In parallel, the conductivity of the solution is measured (Metrohm 712). 

The measurements reveal that all small ions present in solution are free and contribute to the conductivity. This was 
checked at different 7 values by changing the polymer concentration. The measured conductivity is the exact sum of 
the conductivities coming from the Na + ions (whose concentration is known from the amount of added NaOH base) , 
and H + , OH - (known from the pH), while the small contribution of the polymer is negligible. These conductivity 
measurements demonstrate that the distribution of small ions is homogeneous in the solution, and no evidence for 
the Donnan effect and counter-ion condensation is observed. 

IV. RESULTS: COMPARISON OF EXPERIMENTS WITH THEORY 

In this section we present the experimental results for the titration curves of MAA-EA copolymers and compare 
them with the theory of Sec. II. Fig. 1(a) shows the titration curves measured for different polymer concentrations. 
They differ substantially from the titration curves of monomeric methacrylic acid (MAA) . Methacrylic acid in aqueous 
solution has the typical behavior of a weak monomeric acid. The pH increases monotonously with 7, takes the value 
pKa = 4.5 for 7=0.5 and then jumps near 7=1. By contrast, the titration curves of MAA-EA copolymers exhibit 
the following behavior: in the range < 7 < 0.2, the pH increases sharply to a plateau value which remains nearly 
constant up to 7 « 1 where the amount of added base equals that of the monomer (neutralization point). The jump 
of the pH for 7 < 0.2 is associated with a swelling phase transition of the latex particles i.e. a change of conformation 
of the chains. A similar phenomenon has been observed during the titration of pure poly(methacrylic acid) [f23"j . 

Let us mention that it is possible to calculate the swelling transition of the polymer as function of monomer 
concentration, cm. This can be done by minimizing the free energy with respect to c m (in addition to the other 
annealed degrees of freedom discussed above) . Any non-convexity of the free energy signals the existence of a polymer 
precipitate in excess water. 

On Fig. 1 several titration curves are plotted for different concentrations of the same MAA-EA copolymer. As the 
concentration c m increases, the jump in the pH at 7 ~ 1 increases. Note that the deviation between the pH of the 
MAA monomeric acid and the MAA-EA copolymers (c m = 0.1 M in Fig. 1), is large in the plateau region (7 < 1). 
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On the other hand, the deviation is quite small for 7 > 1 when the acid is almost completely dissociated. The value 
of the pH at the plateau depends strongly on the polymer concentration: the larger the c m , the lower the plateau. 
It is surprising that even though MAA-EA copolymers contain carboxylic groups, the plateau value of the pH below 
7 = 1 may be neutral or even greater than 7. 

This difference in behavior between the MAA monomers and MAA-EA copolymers might be associated with a 
complex structure of the MAA-EA copolymer for 7 < 1 . Collapsed microdomains may exist on the chains due to the 
competition between the hydrophobic attraction and Coulomb repulsion. The plateau region suggests the existence 
of such collapsed microdomains along the chains. For 7 1 the chains are in a swollen state and their behavior 
resembles that of monomeric weak acid. 

In Fig. 1(b) we plot for comparison the titration curves as calculated from the theory. The titration equations 
presented in Sec. II. B are solved numerically by an iteration procedure starting from the mean field values as the first 
iteration and including the corrections of the second iteration. Since there are a couple of unknown physical parameters 
(like the second and third virial coefficients, v,w and the monomer size a), we do not try to fit the experimental 
titration curve. Rather, we note that the corresponding theoretical curves look very similar to the experimental ones, 
demonstrating the same type of plateau for low 7 and the same trend with the monomer concentration. 

Figure 2 shows the effect of adding a monovalent salt (NaCl). The overall shape of the titration curves remain 
unchanged. However, the value of the pH at the plateau increases upon the addition of salt, while the pH for 7 > 1 is 
essentially independent on the salt concentration. This indicates that neutralization of the carboxylic groups carried 
by the chains becomes easier as the ionic strength is increased. As the Coulomb interaction is screened by the salt, 
there is less correlation between charged groups along the chain when their distance is larger than the Debye length. 
Therefore, monomers separated by distances larger than the Debye length can dissociate independently from each 
other. For 7 > 1, almost all the charged groups on the chains are dissociated; thus, the salt has no effect on the pH. 
The comparison between experiment and theory is presented in part (a) and (b), where a good agreement can be 
observed. Note the larger deviation from the experimental results for higher salt concentration (cNaCi = 0.1 M). 

In Fig. 3, the plateau value of the pH is plotted on a semi-log plot as function of the ionic strength c\ = CN a + ch + 
c ci + coh- The plateau value of the pH is taken at 7 = 0.5. It is interesting to note that the different points taken 
at various salt and base concentrations collapse on a single curve. This result indicates that the pH is dominated 
by electrostatic effects and not by conformation changes of the chains. When the ionic strength increases, the fixed 
charges carried by the polymer are screened and the energy associated with the dissociation of a carboxylic group 
decreases. 

This behavior can be characterized semi-quantitatively by noting that when 7 ~ 1/2 the pH is much higher than 
the pKa and, therefore, 7 < - - ) — > 1 and A — > 1. Consequently, at the plateau A7 ~ 1/2. To first order in 10 pKA_pH the 
plateau value can be obtained from the simplified free energy, Sec. II. A. 3 and reads: 



P H platcau - P K A + logx 



K/ B + 24^/ 2 c m ^ + 24/^ 



(39) 



Recall that k depends on the total amount of ions in the system through eq. (^|) . The above expression is used for the 

plot of the solid curve of Fig. 3b, where the pH of the plateau is shown as a function of the effective ionic strength 

c G ff defined in eq. ^|. In this regime, the difference between c e g and c\ is a small correction of order l0 pKA ~ pH . 

Equation |3^ exhibits three different regimes depending on the relative importance of the various contributions 

to the logarithmic term. In each regime, ch has a different power law dependence on the total ionic strength c c ff = 

c\ + fA(l — A)c m leading to different slopes in Fig. 3. At low ionic strength the second term dominates and ch oc c~^ 2 

(long dashed line in Fig. 3b). At intermediate ionic strength the last term dominates and ch oc c cff (short dashed 

1/2 

line). Finally, at high ionic strength (beyond the experimentally accessible values) ch oc c cS . 

In Figure 4 the effect of polymer concentration and ionic strength on S(q) is shown. The structure function is 
calculated at the plateau regime were A oil. The dependence on the controlled system parameters: the pH, the salt 
concentration and the degree of neutralization 7 is taken implicitly into account in c e ff. At low values of c e g, the 
structure function exhibits a peak at finite q corresponding to the most favorable density-density fluctuations. As c e ff 
increases, S(q — * 0) increases until the peak disappears. This increase is accompanied by an increase in the osmotic 
compressibility of the solution. 

The effect of varying the polymer concentration c m at a fixed c c ff is depicted in the inset. The peak is stronger at 
lower polymer concentrations indicating that fluctuations become considerably stronger at lower concentrations. This 
effect agrees well with the titration curves shown in Fig. 1. Indeed, at low polymer concentrations the shift in the 
titration curve is stronger than at high concentrations. The peak position shifts to higher wavenumbers (smaller length 



8 



scale) at higher concentration, similar to the correlation length of polymer solutions that decreases with increasing 
polymer concentration. 



V. CONCLUSIONS 



We have presented an experimental study of the titration of weak polyacids by a strong base. Our experimental 
findings performed on solutions of MAA-EA copolymers, are supported by theoretical calculations, which show similar 
trends of the pH variation with the concentration of added base, salt and polyacid. 

Titration experiments are one of the simplest and most useful experimental tools to probe the degree of neutral- 
ization of monomeric and polymeric acids. Since titration curves of weak acids are universal, any deviation from 
this behavior, as observed here, can be the signature of a complex behavior coupling Coulombic and hydrophobic 
interactions. 

The most striking feature is the existence of a plateau of the pH for low degree of neutralization. The plateau in 
the pH is at higher value than the corresponding pH of the monomeric acid indicating that the charges on the chain 
inhibit the dissociation of other charged groups. The pH is not affected strongly by further addition of the base till 
the neutralization point 7 = 1. This is probably due to the existence of collapsed microdomains along the chains 
for 7 < 1. The pH at the plateau decreases as function of the ionic strength. This demonstrates that non-specific 
electrostatic interactions are responsible for the existence of the plateau, since the pH depends mostly on the amount 
of small ions and not on their type (for a fixed degree of neutralization). 

The theory presented above includes one-loop corrections to the mean-field free energy. On the mean-field level, 
there is no difference between the polymeric and monomeric titration curves except for the translational entropy 
of the chains. The one-loop correction couples the chain connectivity with the electrostatics and induces the large 
deviations in the titration curves. We modeled the polyacid as flexible chains using the standard Debye function. Our 
formalism can also be applied to different models of chain elasticity. In particular, for the case of semi-flexible chains 
the same titration equations are obtained with a modified Debye function taking into account the persistence length 
of the chains. 

Finally, it will be interesting to complement this study by scattering experiments where it might be possible to 
resolve the chain microstructures and relate them to the degree of ionization of the chains. 

Acknowledgments: 

IB gratefully acknowledges the support of the Chateaubriand postdoctoral fellowship and the hospitality of the 
Elf-Atochcm research center at Levallois-Perret. DA acknowledges partial support from the US-Israel Binational 
Foundation (BSF) under grant No. 98-00429, and the Israel Science Foundation founded by the Israel Academy of 
Sciences and Humanities — centers of Excellence Program. 



[1] Oosawa, F. Poly electrolytes, Marcel Dekker: New York, 1971. 

[2] De Gennes, P.G. Scaling Concepts in Polymer Physics, Cornell University: Ithaca, 1979. 
[3] Barrat, J.L.; Joanny, J.F. Adv. Chem. Phys. 1996, 94, 1. 
[4] Tanaka,T. Phys. Rev. Lett. 1978, 40, 820. 

[5] Philippova, O.E.; Hourdet, D.; Audebert, R.; Khokhlov, A.R. Macromolecules 1997, 30, 8278. 
[6] Joanny, J.F.; Leibler, L. J. Phys. (France) 1990, 51, 547. 

[7] Kokufuta, E.; Wang, B.; Yoshida, R.; Khokhlov, A.R.; Hirata, M. Macromolecules 1998, 31, 6878. 
[8] Raphael, E.; Joanny, J.F. Europhys. Lett. 1990, 13, 623; Wittmer, J. P. Ph. D thesis (unpublished). 
[9] Hasa, J.; Ilavsky, M. J. Polym. Set. 1975, 13, 263. 
[10] Tirtaatmadja, V.; Tam, K.C.; Jenkins, R.D. Macromolecules 1997, 30, 3271. 

[11] Tam, K.C.; Farmer, M.L.; Jenkins, R.D.; Bassett, D.R. J. Polymer Set.: Part B: Polymer Physics 1998, 36, 2275. 

[12] Katchalsky A.; Michaeli, I. J. Polym. Set. 1955, XV, 69. 

[13] Hasa, J.; Ilavsky, M.; Dusek, K. J. Polym. Sci. 1975, 13, 253. 

[14] Brereton M.G.; Vilgis, T.A. Macromolecules 1990, 23, 2044; Vilgis, T.A.; Borsali, R. Phys. Rev. A. 1991, 43, 6857; 

Brereton, M.G.; Vilgis, T.A. J. Phys. I (France) 1992, 2, 581. 
[15] Napper, D.H. Polymer Stabilization of Colloidal Dispersions, Academic Press: London, 1993. 



9 



[16] Jenkins, R.D.; Lelong, L.M.; Bassett, D.R. in "Hydrophilic Polymers : Performance with Environmental Acceptability" , 

ed. by Glass, J.E., ACS Symposium Series N 248, American Chemical Society 1996, 23, 425. 
[17] Mandel, M. in "Encyclopedia of Polymer Science and Engineering" 1988, 11, 739. 
[18] Strauss, U.P.; Schlesinger, M.S. J. Phys. Chem. 1978, 82, 571. 
[19] Atkins, P.W. Physical Chemistry, Oxford University: Oxford, 1985. 
[20] See, for example, Adv. Polym. Set. 1993, 109; 1993, 110. 
[21] Borukhov, I.; Andelman, D.; Orland, H. Eur. Phys. J. B 1998, 5, 869. 
[22] Auvray, X.; Anthore, R.; Petitpas, C. J. Phys. (France) 1986, 47, 893. 
[23] Heitz, O; Raviso, M.; Frangois, J. Polymer 1999, 40, 1637. 



10 



Fig. 1 pH plotted as a function of the degree of neutralization 7 with no added salt. In (a) the experimental 
titration curves are shown for different polymer concentration: c m = ImM (10~ 4 g/g) full circles; c m = 10 mM (10 -3 
g/g) full triangles; c m = 0.1 M (10~ 2 g/g) full squares. The open circles denote the titration curve of the monomeric 
MAA at c m = 0.1 M (10~ 2 g/g). In (b) the theoretical results are plotted for the same polymer concentrations: 
c m = ImM (dotted curve); c m = 10 mM (short dashes); c m = 0.1 M (long dashes). The solid line corresponds to 
the mean field solution (no correlations and no connectivity) for c m = 0.1 M. It reproduces nicely the conventional 
monomeric titration curve. The polymer parameters we take are: a — 5 A, / = 1/3 and pKa = 4.5. 

Fig. 2 pH plotted as a function of the degree of neutralization 7 for different amount of added salt. The 
experimental values in (a) arc compared with the theoretical ones in (b). The polymer concentration is c m = ImM 
(10~ 4 g/g). The concentration of NaCl in the solution is: CNaCi = (open circles/solid line), 1 mM (full circles/dotted 
curve), 10 mM (full diamonds/short dashes), 0.1 M (full squares/long dashes). The theoretical curves are calculated 
using the same parameters as in Fig. 1. 

Fig. 3 Semi- log plot of the pH at the plateau as a function of the ionic strength. Comparison between experimental 
results (a) and theoretical calculation (b). The experimental measurements are plotted in (a) as function of the total 
amount of small ions c\. The open circles give the pH of salt-free solutions at different monomer concentration 
extracted from Fig. 1. The full symbols give the pH of solutions with increasing amount of added salt: c m = 1 mM 
(10~ 4 g/g) full circles, and c m = 0.1 M (10~ 2 g/g) full triangles. The data collapse indicates that the pH depends 
only on the total ionic strength. The theoretical results are plotted in (b) as a function of c e s = c\ + fA(l — A)c m 
using the same parameters as in Fig. 1 . See text for details of the approximations made in the calculation. 

Fig. 4 The structure function S(q) as function of the wavenumber q. The structure function is calculated at the 
plateau of the titration curve where A ~ 1. The physical parameters used in the calculation are a — 5 A, v = 25 A 3 , 
w = 100 A 6 , TV = 10 4 , f — 1/3 and pKa = 4.5. The monomer concentration is c m = ImM while the effective ionic 
strength is c c e — 0.01 mM (solid line), c e s = ImM (dotted curve) and c c s = 0.1 M (dashed curve). In the inset the 
effective ionic strength is fixed at c c g = ImM while the polymer concentration is varied: c m = ImM (solid line), 
c m = 10 mM (dotted curve) and c m = 0.1 M (dashed curve). 
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Borukhov et al: Fig. 1 
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Borukhov et al: Fig. 3 




Borukhov et al: Fig. 4 
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